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ABSTRACT 

We investigate the nature of molecular gas emission from z ~ 6 quasars via the commonly observed tracer 
of H2, carbon monoxide (CO). We achieve these means by combining non-local thermodynamic equilibrium 
(LTE) radiative transfer calculations with merger-driven models of z ^ 6 quasar formation that arise naturally in 
A-cold dark matter (ACDM) structure formation simulations. Motivated by observational constraints, we con- 
sider four representative 6 quasars formed in the halo mass range ^ 10'^- 10'^ Mq from different merging 
histories. Our main results are as follows. We find that, owing to massive starbursts and funneling of dense gas 
into the nuclear regions of merging galaxies, the CO is highly excited during both the hierarchical buildup of 
the host galaxy and the quasar phase and the CO flux density peaks between J=5-8. The CO morphology of 
z ^ 6 quasars often exhibits multiple CO emission peaks which arise from molecular gas concentrations which 
have not yet fully coalesced. Both of these results are found to be consistent with the sole CO detection at 
z ^ 6, in quasar Jl 148H-5251. Quasars which form at z ^ 6 display a large range of sightline-dependent line 
widths. The sightline dependencies are such that the narrowest line widths are when the rotating molecular gas 
associated with the quasar is viewed face-on (when the Lb is largest), and broadest when the quasar is seen edge 
on (and the Lb is lowest). Thus, we find that for all models selection effects exist such that quasars selected 
for optical luminosity are preferentially seen to be face-on which may result in CO detections of optically lu- 
minous quasars at z ^ 6 having line widths narrower than the median. The mean sightline-averaged line width 
is found to be reflective of the circular velocity of the host halo, and thus scales with halo mass. For example, 
the mean line width for the ^ 10'^ Mq halo is cr ^ 300 km s"', while the median for the ~ lO'^* Mq quasar 
host is a ^ 650 km s"'. Depending on the host halo mass, approximately 2-10% of sightlines in our modeled 
quasars are found to have narrow line widths compatible with observations of J1148H-5251. When considering 
the aforementioned selection effects, these percentages increase to 10-25% for quasars selected for optical lu- 
minosity. When accounting for both temporal evolution of CO line widths in galaxies, as well as the redshift 
evolution of halo circular velocities, these models can self-consistently account for the observed line widths 
of both submillimeter galaxies and quasars at z 2. Finally, we find that the dynamical mass derived from 
the mean sightline-averaged line widths provide a good estimate of the total mass, and allows for a massive 
molecular reservoir, supermassive black hole, and stellar bulge, consistent with the local MaH-A^buige relation. 

Subject headings: cosmology: theory — cosmology: early universe — galaxies: formation — galaxies: active 
— galaxies: high-redshift — galaxies: ISM — galaxies: individual (SDSS J1148H-5251) 



1. INTRODUCTION 

The discovery of extremely luminous quasars at z >6 via 
novel color selection techniques demonstrates that massive 
galaxies and supermassive black holes formed very early in 
the Universe (Fan et al. 2002, 2003, 2004). Observations of 
the dusty and molecular interstellar medium (ISM) in z ^ 6 
galaxies can serve as a unique probe into the star formation 
process in the first collapsed objects, and help quantify the re- 
lationship between star formation and black hole growth when 
the Universe was less than a billion years old (e.g. Wang et 
al. 2007; for a recent review, see Solomon & Vanden Bout 
2005). Serving as a proxy for observationally elusive molec- 
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ular hydrogen (H2), rotational transitions in tracer molecules 
such as '^CO (hereafter, CO), HCN and HCO+ can provide 
diagnostics for the physical conditions in the star-forming gi- 
ant molecular clouds (GMCs) of high redshift galaxies (e.g. 
Bertoldi et al. 2003a,b; CarilH et al. 2005; Riechers et al. 
2006a; Walter et al. 2003). 

The highest redshift quasar that has been found, SDSS 
J1148H-5251 (hereafter J1148H-5251) at z=6.42 (Fan et al. 
2003) is an extremely bright object with bolometric luminos- 
ity lO'"* Lq, and is thought to be powered by accretion 
onto a supermassive black hole of mass ^10^ Mq (Willott, 
McLure, & Jarvis 2003). Bertoldi et al. (2003b) measured a 
far infrared luminosity of 1.3 x 10'^ L©, which, if powered 
solely by starburst-heated dust, corresponds to an exceptional 
star formation rate of ^ 3000 Moyr"'. 

Pioneering CO observations of Jl 148H-525 1 have revealed a 
great deal concerning the molecular ISM in the host galaxy of 
this z ^ 6 quasar. Through multi-line observations and large 
velocity gradient (LVG) radiative transfer modeling, Bertoldi 
et al. (2003a) found that the CO flux density peaks at the J=6 
level of CO, indicative of the warm and dense conditions char- 
acteristic of vigorous star formation. A measured molecular 
gas mass of Mh^ « 10'" Mq shows that J1148H-5251 plays 
host to a large reservoir of molecular gas (Bertoldi et al. 
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2003a; Walter et al. 2003). Subsequent high resolution ob- 
servations with the Very Large Array (VLA) by Walter et 
al. (2004) discovered that the CO emission in this galaxy 
is extended on scales of 2.5 kpc and is resolvable into two 
emission peaks separated by 1.7 kpc, with each peak trac- 
ing ~ 5 X lO' Mq of molecular gas. These observations sug- 
gest that ni48H-5251 is a merger product (Walter et al. 2004; 
Solomon & Vanden Bout 2005). 

CO observations of Jl 148H-5251 have presented challenges 
for models of galaxy formation as well. For example, dynami- 
cal mass estimates from observed CO line widths are not able 
to account for the presence of a ^ lO'^ M© stellar bulge as 
would be predicted by the present-day MsH-Mbuige relation, 
suggesting that the central supermassive black hole could 
have grown in part before the host galaxy (Walter et al. 2004). 
In contrast, the presence of heavy elements (Barth et al. 
2003), and significant CO emission (Bertoldi et al. 2003a; 
Walter et al. 2003, 2004) imply that the ISM has been sig- 
nificantly enriched with metals from early and abundant star 
formation. Recent theoretical arguments have additionally 
proposed that a relation between black hole mass and stel- 
lar bulge mass is a natural consequence of AGN feedback in 
galaxies (Di Matteo et al. 2005, 2007; Hopkins et al. 2007a; 
Sijacki et al. 2007), and that this relation shows only weak 
(~ 0.3-0.5 dex) evolution in galaxies from redshifts z=0-6 
(Robertson et al. 2006a,c; Hopkins et al. 2007a). 

Numerical simulations can offer complementary informa- 
tion to observations of z ^ 6 quasars by providing a frame- 
work for the formation and evolution of z ^ 6 quasars and 
the relationship between the star forming ISM in these galax- 
ies and observed CO emission. Calculations by Springel, Di 
Matteo & Hernquist (2005a,b) have found that galaxy merg- 
ers serve as a viable precursor for quasar formation. Strong 
gaseous inflows driven by tidal torques on the gas (e.g. Barnes 
& Hernquist 1991, 1996) can fuel nuclear starbursts (e.g. Mi- 
hos & Hernquist 1994, 1996) and feed the growth of central, 
supermassive black holes (Di Matteo et al. 2005); subse- 
quent feedback from the AGN can lift the veil of obscuring 
gas and dust and, along numerous sight lines, reveal an opti- 
cally bright quasar (Hopkins et al. 2005a,b). 

More recently, Li et al. (2007) have proposed a merger 
driven model for quasar formation at z ^ 6 which fits nat- 
urally into a ACDM framework. By performing numeri- 
cal simulations which simultaneously account for black hole 
growth, star formation, quasar activity, and host spheroid 
formation, these authors found that galaxy mergers in early 
~ 10^^ M0 halos can result in the formation of bright quasars 
at z ~ 6. The quasars in these simulations exhibit many prop- 
erties similar to the most luminous quasars at z ^ 6, including 
both observables, such as the rest-frame £-band luminosity, 
and inferred characteristics (e.g. the central black hole mass 
and bolometric luminosity). These simulations can thus serve 
as a laboratory for investigating the properties of the interstel- 
lar medium in the first quasars, as well as the relation between 
CO emission and star formation at early epochs. 

In order to quantitatively couple these models with CO ob- 
servations, molecular line radiative transfer calculations are 
necessary. Recent works by Narayanan et al. (2006a,b) have 
developed a methodology for simulating molecular line trans- 
fer on galaxy-wide scales. Here, we aim to investigate the 
plausibility of z ~ 6 quasar formation in massive ~ lO'^-lO^^ 
Mq halos in a hierarchical structure formation scenario by 
modeling the observed CO emission from these high redshift 
sources. We achieve these means by coupling the non-LTE 



radiative transfer codes of Narayanan et al. (2006a,b) with 
the hierarchical z ^ 6 quasar formation models of Li et al. 
(2007). We make quantitative predictions for z ~ 6 quasars 
which form in halos with masses ranging from ^ 10'^ - 10^^ 
Mq, and provide interpretation for existing and future CO ob- 
servations of z >6 quasars. 

In § 2 we describe the numerical models employed. In § 3, 
we present the modeled CO excitation characteristics and lu- 
minosities. In § 4, we discuss the morphology of the CO gas. 
In § 5 we describe the derived emission line profiles and com- 
pare our models to observations of the only detected CO mea- 
surement in a z ^ 6 quasar, J1148H-5251. In § 6, we use the 
model CO lines to investigate the usage of CO observations 
as dynamical mass indicators in the highest redshift quasars. 
We conclude with a discussion comparing the CO emission 
properties of our simulated quasar with other high-z galaxy 
populations in §7 and summarize in § 8. Throughout this pa- 
per we assume a cosmology with h=OJ, Oa=0.7, ^Im=0.3. 

2. NUMERICAL METHODS 

In order to capture the physics of early Universe quasar for- 
mation, simulations must have the dynamic range to faithfully 
track the evolution of the most massive halos within which 
these quasars reside (e.g. Haiman & Loeb 2001), as well as 
foUow the stars, dark matter, ISM and black holes in the pro- 
genitors of the quasar host galaxy. We have performed multi- 
scale calculations which include cosmological dark matter 
simulations in a volume of 3 Gpc^ to identify the most mas- 
sive halos, and subsequent hydrodynamic galaxy merger com- 
putations within these halos at higher resolution. The cos- 
mological and galaxy merger simulations were implemented 
using the parallel, A'-body/Smoothed Particle Hydrodynam- 
ics (SPH) code GADGET-2 (Springel 2005). We then appHed 
the non-LTE radiative transfer code, Turtlebeach (Narayanan 
et al. 2006a,b), to the outputs of the hydrodynamic galaxy 
merger simulations to investigate the emission properties of 
the CO molecular gas in the resultant quasar The quasar 
formation simulations and radiative transfer algorithms are 
described in detail in Li et al. (2007) and Narayanan et al. 
(2006a,b), respectively, and we refer the reader to those works 
for more detail; here, we briefly summarize, and focus on the 
aspects of the calculations directly relevant to this study. 

2.1. Cosmological Simulations 

Observational evidence suggests that quasars at z ~ 6 Ukely 
form in very massive (~ 10^^ Mq) halos. For example, 
quasars at z ^ 6 have an extremely low comoving space den- 
sity of n « 10"'^ Mpc~^ (e.g. Fan et al. 2003), compara- 
ble to the rarity of massive ~ lO'^ Mq halos at z ^ 6. Fur- 
ther evidence comes from black hole mass-halo mass corre- 
lations, and theoretical arguments relating quasar luminosity 
and halo mass (Lidz et al. 2006). In addition, recent numeri- 
cal simulations by Pelupessy, Di Matteo & Ciardi (2007) have 
suggested that the buildup of supermassive ~ 10^ Mq black 
holes by z ~ 6 may require galaxy mergers in extremely rare 
> 10^^ Mq halos. The strongest evidence for high redshift 
quasar formation in massive halos comes from recent SDSS 
clustering measurements by Shen et al. (2007) which indi- 
cate that the minimum mass for high redshift quasar hosts is 
~6-9x10'^Mq. While these arguments are suggestive, 
they are not conclusive. Thus, one major goal of this work is 
to further investigate the hypothesis that z 6 quasars live in 
high mass halos and constrain the range of host halo masses 
for these sources by comparing their simulated CO emission 
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to observations. To this end, we have simulated 4 quasars 
which form hierarchically through numerous mergers in mas- 
sive halos ranging in mass from ~ 10'^- 10^^ M©. 

In order to simulate the large subvolume of the Universe 
necessary to track the evolution of massive halos and forma- 
tion of z ~ 6 quasars, as well as achieve a reasonable mass 
resolution at the same time, we performed the structure for- 
mation simulation with a multi-grid procedure similar to Gao 
et al. (2005). We first ran a uniform resolution simulation of a 
3 Gpc^ volume, with an effective mass and spatial resolution 
of mdm ^ 1.3 X 10^-^ M0 and e ~ 125 h~'kpc (comoving soft- 
ening length), respectively, with initial conditions generated 
by CMBFAST (Seljak & Zaldarriaga 1996). We assumed a 
(78 of 0.9 (though see Li et al. [2007] for a discussion of the 
impUcations of other choices for the cosmological parame- 
ters). 

We then used a friends-of-friends group finder in order to 
seek out candidate massive halos within which early quasars 
form. Anticipating that these rare quasars will be progenitors 
of massive objects today, we identified the most massive halos 
at z=0, and then resimulated the evolution of these objects 
and their immediate environment at a much higher mass and 
force resolution assuming an initial redshift of z=69. The final 
effective resolutions of the halo evolution simulations were 
OTdm ~ 2.8 X 10*^ Mq and e ~ 5 h'^kpc . 

The merger tree of the halos was extracted from the simula- 
tions in order to provide information concerning the masses of 
the largest contributing progenitors to the halo mass and allow 
us to reconstruct the hierarchical buildup of the quasar host 
galaxy. We considered groups that contribute at least 10% of 
the halo mass as progenitors in the merger history resulting in 
e.g. seven progenitors for the most massive halo simulated. 
Table 1 lists the properties of all the resultant quasars. 

The discussion throughout this paper will largely focus on 
quasars Q1-Q3 which formed hierarchically in the cosmolog- 
ical simulations. Quasar Q4 is a binary coplanar merger spe- 
cially simulated for studying the dependence of the CO emis- 
sion properties on merger history, and will be discussed only 
in §5.2. 

2.2. Galaxy Merger Simulations 

To derive the physical properties of the z ^ 6 quasar host 
galaxies, we resimulated their merger trees hydrodynamically 
using GADGET-2. GADGET-2 utilizes a fully conservative 
SPH formulation which allows for an accurate handUng of 
discontinuities (Springel & Hemquist 2002). The code ac- 
counts for radiative cooling of the gas (Katz et al. 1996; Dave 
et al. 1999), and a multiphase description of the ISM that in- 
cludes cold clouds in pressure equilibrium with a hot, diffuse 
gas (e.g. McKee & Ostriker 1977; see also Springel & Hern- 
quist 2003a). Star formation is constrained by observations of 
local galaxies, and follows the Kennicutt-Schmidt laws (Ken- 
nicutt, 1998; Schmidt 1959; Springel & Hernquist 2003b). 
The progenitor galaxies had dark matter halos initialized to 
follow a Hernquist (1990) profile, and the virial properties are 
scaled to be appropriate for cosmological redshifts (Robertson 
et al. 2006a). 

Black holes in the simulations are realized through sink par- 
ticles that accrete gas following a Bondi-Hoyle-Lyttleton pa- 
rameterization (Bondi & Hoyle 1944; Bondi 1952; Hoyle & 
Lyttleton 1939). To model feedback from central black holes, 
we assume that 0.5% of the accreted mass energy is rein- 
jected into the ISM as thermal energy (Di Matteo et al. 2005; 
Springel, Di Matteo & Hemquist 2005a,b). This formula- 



tion for AGN feedback in galactic scale simulations has been 
shown to successfully reproduce X-ray emission patterns in 
galaxy mergers (Cox et al. 2006a), observed quasar luminos- 
ity functions and lifetimes (Hopkins et al. 2005a-d, 2006a,c,d, 
2007b), the Seyfert galaxy luminosity function (Hopkins & 
Hernquist 2006), the Mbh-o-v relation (Di Matteo et al. 2005; 
Hopkins et al. 2007a; Robertson et al. 2006a), the bimodal 
galaxy color distribution (Springel et al. 2005b, Hopkins et 
al. 2006b,e, 2007c), characteristic CO emission patterns in 
ultraluminous infrared galaxies (ULIRGs; Narayanan et al. 
2006a), infrared colors of ULIRGs and z ~ 2 Submillimeter 
Galaxies with embedded AGN (Chakrabarti et al. 2007a,b), 
and the kinematic structure of merger remnants (Cox et al. 
2006b). 

We assume that the black holes in these simulations formed 
from the first stars (e.g. Abel, Bryan & Norman 2002; Bromm 
& Larson 2004; Yoshida et al. 2006), and that the seed 
masses for the black holes are 200 M0 at z=30. Before the 
first progenitor galaxies entered the merger tree, their black 
holes grew at the Eddington Umit, resulting in seed black hole 
masses of ~ 10^ Mq at the time of the first merger. The black 
holes in the simulations are assumed to merge when their sep- 
aration is less than a gravitational softening length (30 h"'pc). 

We model evolution of the CO emission in the quasars 
from z w8 to z «5, noting that the peak of the quasar activ- 
ity is roughly 7> z ^5.5. The black hole luminosity outshines 
the stellar luminosity for a large range of redshifts, and it is 
this time that we refer to as the 'quasar phase' for any given 
galaxy. 

2.3. Evolution of Model Quasars 

To aid in the discussion in the remainder of this work, 
we qualitatively describe the evolution of the most massive 
(^ 10'^ M0) quasar host, though the results are general for 
all models considered in Table 1. The quasar host galaxy 
builds hierarchically, through seven major mergers between 
z= 14.4-8.5. Strong gravitational torques on the gas drive 
massive gaseous inflows, causing heavy accretion onto the 
central black hole(s), and triggering intense starbursts that 
typically form stars at a rate between ^ 10^- 10"* MQyr"'. 
The black holes accrete heavily as gas is funneled in toward 
the nuclear regions. Feedback from the most massive central 
black hole then drives a powerful wind, creating numerous 
lines of sight along which the central quasar is no longer ob- 
scured, and the black hole luminosity outshines the stellar lu- 
minosity. The central supermassive black hole can be visible 
as an optically bright quasar (Lboi > IO^^Lq) for 50 Myr 
(though less for lower mass quasar host galaxies; Table I). 
The powerful quasar wind quenches the starburst, and self- 
regulates the black hole growth. In the post quasar phase, the 
luminosity of the galaxy subsides, and it eventually evolves 
into a cD-like galaxy. More details concerning the evolution 
of the models presented here are discussed in Li et al. (2007). 

2.4. Molecular Line Radiative Transfer 

The CO emission properties were calculated using Turtle- 
beach, a 3-dimensional non-LTE radiative transfer code, 
based on an expanded version of the Bernes (1979) algorithm 
(Narayanan et al. 2006a,b). Our improvements focus on in- 
cluding a mass spectrum of GMCs in a subgrid manner, which 
allows us to more accurately model the molecular line radia- 
tion from the dense cores of molecular clouds, as well as the 
diffuse outer layers. 
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TABLE 1 

Quasar Models 



Name" z" Mvi,' Vvir"' Mbh^ iboi(peak)f TQso(Lboi > 10^^ Lq)^ 

[lO'^M©] [kms"'] [10' Mq] [10" L©] [Myr] 



Ql 6.5 7.7 626.4 2.0 2.0 55.5 

Q2 5.6 3.8 467.9 1.0 5.0 22.3 

Q3 7.2 1.5 340.7 0.2 1.5 4.0 

Q4 5.5 8.1 593.2 0.4 1.1 21.9 



^Name of quasar model. 

^Redshift of the quasar at peak accretion activity. 

Virial mass of the quasar host halo, assuming overdensity A = 200. 

Virial velocity of the quasar host halo, assuming overdensity A = 200. 
^Black hole mass of the quasar. 
^Peak bolometric luminosity of the quasar. 
^Quasar hfetime for Lboi > 10 Lq. 



We build the emergent spectrum by integrating the equation 
of radiative transfer along various Unes of sight: 

r 

ro 

where I^y is the frequency-dependent intensity, 5,y is the source 
function, r is the physical depth along the line of sight, and r 
is the optical depth. The source function is dependent on the 
CO level populations which are assumed to be in statistical 
equilibrium, but not in LTE. This aspect of the calculation is 
important as the assumption of LTE often breaks down when 
considering the propagation of submilUmeter-wave radiation 
through the molecular ISM. 

We consider collisions with H2, radiative excitation, and 
stimulated and spontaneous emission in determining the CO 
level populations. The radiative transfer is handled in a Monte 
Carlo manner, in which photon 'packets' are emitted isotropi- 
caUy over47r steradians. The frequencies are randomly drawn 
from the local line profile function which includes effects 
from the kinetic temperature and a microturbulent velocity 
field which we assume to have a constant value of 0.8 km s~' . 
The strongly density-dependent collisional excitation rates are 
modeled by deriving the density distribution in a given cell 
from a mass spectrum of giant molecular clouds. In this for- 
mulation, we model the density distributions in the clouds as 
singular isothermal spheres with power-law index p=2 (e.g. 
Walker, Adams & Lada 1990), where the radius of each cloud 
is determined by the Galactic GMC mass-radius relation (e.g. 
Solomon et al. 1987; Rosolowsky 2005, 2006). Our results 
for higher lying CO transitions are particularly improved us- 
ing this subgrid non-LTE approach. The molecular mass frac- 
tion of cold gas is assumed to be half, as motivated by local 
volume surveys (e.g. Keres, Yun & Young, 2003). This re- 
sulted in molecular gas masses of 1 - 3 x IO^^Mqui the model 
quasars. 

Our methodology is an iterative one. A solution for the CO 
level populations is first guessed and provides the initial ra- 
diation field. We then calculate the mean intensity field in a 
Monte Carlo manner. The contributions of the raiation field 
and collisional excitation rates determine the updated level 
populations via rate equations and the process is repeated un- 
til convergence. 

We have benchmarked our radiative transfer codes against 
the published non-LTE line transfer tests of van Zadelhoff et 
al. (2002), and present the results in Narayanan et al. (2006b). 
Our approach in combining non-LTE radiative transfer calcu- 
lations with galaxy-scale hydrodynamic simulations has pre- 



dicted CO morphologies, spatial extents, excitation levels, 
line widths and intensities that are in good agreement with ob- 
servations of local starbursts and galaxy mergers (Narayanan 
et al. 2006a; D. Narayanan et al, in prep.). 

For the models presented here, 1x10^ model photons were 
emitted per iteration. The mass spectrum of GMCs in each 
sub-grid cell considered clouds with a lower cutoff at 1 x 10"* 
M0 and an upper cutoff at 1 x 10'' Mq, consistent with ob- 
served mass ranges of GMCs (e.g. Blitz et al. 2006). We cal- 
culated CO transitions across 10 levels at a time and assumed 
a uniform Galactic CO abundance of CO = 1.5 x 10"^ x H2 
(Lee, Bettens & Herbst 1996). We utilized collisional rate 
coefficients from the Leiden Atomic and Molecular Database 
(Schoier et al. 2005). 

Throughout this work, we make predictions for CO tran- 
sitions ranging from J=l-0 through J=10-9. One potential 
caveat regarding the lower lying (e.g. J=l-0) transitions is 
that the low density molecular gas in z 6 quasars may be 
in a diffuse intracloud medium as is thought to be the case 
for local ULIRGs, rather than bound in GMCs (Downes & 
Solomon 1998). In our subgrid modeling of GMCs, we as- 
sume implicitly that the diffuse molecular gas is bound in the 
outer envelopes of GMCs which are characterized as singular 
isothermal spheres. The assumption of higher density gas re- 
siding in bound cloud cores is likely to be more robust. Recent 
studies show similarities between LIR/HCN ratios in both lo- 
cal GMCs and a large sample of ULIRGs. These similarities 
suggest that the high density gas even in the nuclei of galaxy 
mergers is locked in cloud cores (Wu et al. 2005), similar to 
the simulations presented here. In this sense, our modeled en- 
vironment of higher density gas is more reasonable, and pre- 
dictions concerning higher-lying transitions of CO (e.g. J>3) 
more robust. 

3. EXCITATION AND LUMINOSITY OF CO 

The CO spectral energy distribution (SED) is the CO flux 
density emitted per rotational J state and thus provides a di- 
rect measure of the excitation characteristics of molecular gas 
(e.g. WeiB et al. 2005a). The shape and peak of the CO SED 
describe the relative number of molecules in a given rotational 
J state, and serve as observable diagnostics of the underlying 
temperature and density of the emitting molecular gas. Here, 
we describe in detail the excitation characteristics of the most 
massive quasar host (Ql, M ^ lO*-' Mq), and expUcitly note 
when the lower mass models exhibit different properties. In 
Figure 1, we show the redshift evolution of the peak in the 
CO SED. Moreover, in Figure 2, we show representative CO 
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Fig. 1. — Upper rotational state (Jupp^r) of CO SED peak as a function of 
redshift for most massive quasar host, model Ql (see Table 1). During the 
early, massive starburst, most of the CO is highly excited. During the height 
of the quasar phase, the peak in the CO excitation in our simulated galaxy 
ranges from J=8 in the beginning of the quasar phase to J=5 near the end. 
This is broadly consistent with observations of excited molecular gas in the 
sole CO detection at z ~ 6 (Bertoldi et al. 2003a). As the starburst subsides 
in the post-quasar phase, lower temperatures and densities drive the peak in 
the CO SED down to J ~ 3. 

SEDs for quasar Ql at three redshifts during the quasar phase. 
The CO SEDs are averaged over three orthogonal viewing an- 
gles, though the SEDs derived from each individual angle are 
nearly identical. 

Strong gravitational torques exerted on the gas by the multi- 
ple mergers drive massive amounts of cold gas into the central 
kiloparsec, giving rise to densities as high as ^ 10^ cm"-' in 
GMC cores at the beginning of the quasar phase (z ^ 7). The 
combination of these dense conditions and heating associated 
with the continued starburst cause the CO molecular gas to 
become highly excited. During this time, the peak of the CO 
SED rises to J=8. Molecular gas in the central ~ kiloparsec 
dominates the high-lying excitation. To place this in the con- 
text of active star forming regions in the local Universe, WeiB, 
Walter & Scoville (2005b) found the CO SED to peak at J=6 
in the nuclear region of starburst galaxy M82. 

As energy input from the quasar begins to quench the star- 
burst and dispel gas from the central regions, the star forma- 
tion rate drops to ~ 10^ M©yr~^ Consequendy, the peak ex- 
citation in the molecular gas rapidly drops to values more sim- 
ilar to local starbursts. During most of the quasar phase, the 
CO flux density in the simulated galaxy peaks at J=6. This 
is consistent with the multi-hne observations and CO SED 
derivations by Bertoldi et al. (2003a), who find the CO flux 
density in J1148H-5251 to peak at J=6. In a merger-driven 
model for high-z quasars, the interplay between massive star- 
bursts and feedback from central black holes is important 
in determining the observed CO excitation characteristics in 
galaxies Uke J1148H-5251. 

The peak of the CO flux density remains roughly constant 
through the height of the quasar phase. As feedback from 
the central black holes further extinguishes the nuclear star- 
burst, fewer molecules are highly excited, and consequently, 
the relative flux density from higher J levels begins to drop. 
In Figure 2, the slope of the CO SED at levels higher than the 
turnover point is seen to become steeper as the quasar evolves. 
When the accretion onto the central supermassive black hole 
subsides in the post-quasar phase (z <6), the star formation 
rate (SFR) drops to <50 Mq yr"^ The bulk of the molecular 
gas in this late stage of the galaxy's evolution is only moder- 
ately excited, and the peak in the CO SED declines to J«3-4. 

The trends discussed above are similar in models Q2 and 



Q3 (Table 1), though the overall normalization is slightly dif- 
ferent. The CO SED of both the intermediate mass and low 
mass models (Q2 and Q3) peaks at J=6 at the beginning of 
the quasar phase, and settles at J=4 as the starburst subsides, 
phase. The lower excitation values in the lowest mass model 
owe to overall lower densities and star formation rates. For 
example, during the quasar phase, the SFR from model Q3 is 
~30 Moyr"^ 

Another way to view the excitation characteristics of the 
molecular gas is through the velocity-integrated CO luminos- 
ity (in units of K-km s~\ where the K is the Ray leigh- Jeans 
temperature). In Figure 3, we show the normalized veloc- 
ity integrated CO (J=l-0, J=3-2, J=6-5 and J=9-8) intensity 
as a function of redshift. We additionally plot the SFR, stel- 
lar, and black hole luminosity. The CO luminosity across all 
transitions decreases as the merger activity progresses, and 
the starburst reaches its peak. This is especially true of the 
high-J states which are typically excited by collisions in the 
starburst-heated gas. All CO transitions peak in integrated in- 
tensity early on, when the starburst has not consumed most 
of the available star-forming gas, and collisions help to sus- 
tain molecular excitation. As the starburst fades owing to a 
diminishing fuel supply, the intensity from the high lying CO 
transitions (e.g. 6-5, 9-8) faUs off rapidly while the lower-J 
transitions experience a more moderate decline. In part, this 
owes to the fact that the lowering of gas temperatures and den- 
sities does not heavily affect the molecular excitation at J=l. 
De-excitation of warm, dense star forming gas additionally 
contributes to populating the lower J states. While the CO lu- 
minosity is only about half of its maximum value during the 
quasar phase, the bolometric luminosity of the galaxy peaks 
here as the central quasar becomes visible (Figure 3, bottom 
panel). 

4. MOLECULAR GAS MORPHOLOGY 

In this section, we discuss the CO morphology of the model 
quasar host galaxies through their evolution. The discussion 
is again focused on the most massive simulation (Ql), though 
the results are generic for each of the models studied here. 
In Figure 4, we show the evolution of the central 2 kpc of 
the CO (J=l-0) emission in the most massive halo, Ql, dur- 
ing the hierarchical buildup of the host galaxy and quasar 
phase. Individual concentrations in the molecular gas density 
which have not fully coalesced appear during the buildup of 
the quasar host galaxy, and through parts of the quasar phase 
(e.g. z=6.73. Figure 4), giving rise to multiple CO (J=l-0) 
emission peaks. Near the end of the quasar's lifetime, the 
molecular gas settles into a nuclear disk, with the densest gas 
in the central ~ 500 pc. 

The morphological features of our simulated quasar 
agree reasonably well with observations. Observations of 
Jl 148H-5251 have revealed two CO (J=3-2) emission peaks in 
the central 2 kpc (Walter et al. 2004), similar to the multi- 
ple surface CO surface brightness peaks seen at many points 
in our models (e.g. Figure 4). This suggests that the ob- 
served multiple surface brightness peaks in the CO morphol- 
ogy of Jl 148H-525 1 may owe to separated peaks of high den- 
sity emission in the nucleus that have not yet coalesced. To 
further illustrate this point, in Figure 5, we have plotted the 
CO (J=3-2) emission contours at z=6.73 over three orthogo- 
nal viewing angles, and six random ones. Most viewing an- 
gles exhibit two surface brightness peaks, suggesting that a 
merger origin for the formation of Jl 148H-5251 is viable. 

Within the constraints of our numerical simulations, multi- 
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Fig. 2. — The CO SED (solid line) at three points during the peak of the quasar phase of our most massive simulated galaxy (M ~ lO'^ Mq; model Ql). 
Massive starbursts and dense conditions in the central kiloparsec cause the CO flux density to peak at J=8 prior to and at the beginning of the quasar phase (not 
shown). During the bulk of the quasar phase, the CO SED peak drops to J=6 owing to quenching of the starburst by black hole feedback, and remains roughly 
constant for the remainder of the quasar's life. As the central AGN quenches the nuclear starburst, less gas is in highly excited (J ^6) states, and the slope on the 
blue end of the peak becomes much steeper. The dotted line is the CO SED from just the central kpc, and the dashed line is the contribution from radius=l-4 
kpc). The upper level of each transition is listed on the top axis. 

bursts can cause transition-dependent CO morphologies. In 
Figure 6, we show the CO (J=l-0), (J=3-2), and (J=6-5) emis- 
sion contours for the central 2 kpc of the most massive host 
during the quasar phase. The multiple emission peaks which 
owe to merging clumps of dense gas (z=6.73 in this plot) ap- 
pear in all transitions. When the cold gas has coalesced in 
the nucleus, emission from the lowest excitation gas traced 
by CO (J=l-0) exhibits discrete pockets of emission whereas 
higher-lying emission originating from denser gas is centrally 
concentrated. As seen in Figure 2, the warm and dense condi- 
tions in the starbursting nucleus of the quasar ensure that the 
CO remains excited, with the peak flux density at J=5-8. The 
CO gas in the nucleus is typically excited out of the J=l level, 
resulting in discrete pockets of emission in the circumnuclear 
molecular gas. The higher-lying emission, which originates 
in the warmer, denser cores of GMCs, is more prevalent in 
the central 500 pc of the quasar, resulting in a smoother dis- 
tribution. 

Owing to the violent nature of the quasar's formation pro- 
cess, the spatial extent of the CO emission has a large dynamic 
range, and varies from ~ 2 kpc to ~ 300 pc (half-light radius) 
throughout the evolution of the host galaxy. The more ex- 
tended morphologies are representative of times when gas is 
falling into the central potential, and in the post quasar phase 
when winds have produced an extended morphology. More 
compact CO emission is seen primarily when the cold gas has 
completely coalesced, and the quasar is most active. 

The starburst is centrally concentrated in the central ^ 250 
pc of the quasar host galaxy during the majority of the sim- 
ulation presented here. Conversely, the CO morphology is 
not always so compact. The CO emission through J=10 is 
extended during the buildup of the host galaxy (prior to the 
quasar phase) owing to merging gas clumps. This implies 
that CO may not always serve as an adequate tracer of the 
starburst during the major merger phase of the quasar's for- 
mation. During the active quasar phase, when most of the gas 
has fallen into the nucleus, emission from higher CO transi- 
tions (e.g. J>6) tends to become compact, and faithfully trace 
the active starburst. 

5. CO EMISSION LINES 
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Fig. 3. — Top: Velocity-integrated CO Intensity (Kkm s~') as a function of 
redshift for various CO transitions for quasar Ql. Each transition is normal- 
ized. Middle: Star formation rate and Bottom Quasar bolometric luminosity 
as a function of redshift. The bolometric luminosity is broken up into contri- 
bution by stars and black holes. The bright CO emission in high J transitions 
dming the quasar phase is representative of massive star formation, though 
the bulk of the spheroid formed during the hierarchical buildup of the quasar 
host galaxy, prior to the active quasar phase. 

pie density peaks in the cold gas appear to be the most plausi- 
ble explanation for the observed morphology of Jl 148-1-5251. 
Multiple emission peaks in the CO morphology of galaxy 
mergers have also been noted to arise from large entrainments 
of molecular clouds in AGN-driven winds (Narayanan et al. 
2006a). In the current simulations, however, the character- 
istic outflow morphologies of Narayanan et al. are not seen 
during the active quasar phase. Emission from the nuclei of 
progenitor galaxies can additionally cause multiple CO sur- 
face brightness peaks. However, by the time the simulated 
host galaxy reaches the height of the quasar phase, the most 
massive nuclei have all merged into the central potential (Li 
et al. 2007). 

The excitation characteristics in the vicinity of massive star- 



The simulated properties of the CO emission lines show 
variations among the models presented in Table 1. In this 
section we discuss the CO emission Unes in terms of the CO 



The Nature of CO Emission in z ~ 6 Quasars 



7 




z=6.73 t:' 




z=6.54 



9=D 



z=6.03 



z=5.88 



9=D U_ 



z=5.74 « 

l|)=D 



I 



Te (K kin s"0 



2L79 



4358 



6538 



8717 



Fig. 4. — Evolution of CO (J=l-0) emission contours from z=8.2-5.7 in most massive halo, model Ql (M ^ lO'^ Mq). Multiple emission peaks are visible as 
cold gas falls in toward the nucleus, similar to observations of J1148+5251. As the starburst and quasar activity subsides, the CO intensity fades. Each panel is 
four kpc on a side. The emission is in terms of velocity-integrated intensity (K-km s"'), and the scale is at the bottom. The viewing angle of each panel is listed 
in the top right corner, and is always = 0,4> = for this figure. 



(J=6-5) lines as the CO SED peaks near J=6 for most of the 
quasar phase, and thus this transition best traces the properties 
of the bulk of the molecular gas. Unless otherwise specified, 
the nature of the CO line profile as described in the remainder 
of this section is not seen to vary significantly with observed 
transition. 

5.1. General Nature of Modeled Line Profiles 

In Figure 7, we show a sample 3 random CO (J=6-5) emis- 
sion lines from the simulated quasars Ql, Q2 and Q3 at the 
peak of their respective quasar phases. The first noteworthy 
point regarding the CO emission lines from the 6 quasar 
models is that they are characteristically featured, with a sig- 
nificant amount of substructure. Along many lines of sight, 
several smaller emission spikes (with a 50- 100 km s"') 
originating in dense CO clumps in the central 2 kpc sit super- 
posed on the broader emission line. These emitting clumps 
are centrally concentrated in the quasar, and similar emission 
features are seen in higher spatial resolution spectra. 

Another clear feature of Figure 7 is the apparent trend of 
narrowing line width with decreasing halo circular velocity. 
The CO linewidths are reflective of the circular velocity, and 
are thus strongly dependent on the mass of the quasar host 
galaxy. To further illustrate this, in Figure 8, we show the 
mean sightline-averaged CO (J=6-5) line width during the hi- 
erarchical buildup and quasar phase of the quasars formed 
in the halos of Q1-Q3. We additionally show the range of 
sightline-dependent observed line width values. The multi- 
ple mergers involved in the formation of the simulated 6 
quasars give rise to large velocity dispersions along a num- 
ber of sightlines during the buildup of the quasar host galax- 
ies. During this time, much of the molecular gas is not virial- 
ized, and thus the typical line widths exhibited represent about 
twice the expected circular velocity at the spatial extents of 
the molecular gas. In the quasar phase, as the gas virializes 



into a molecular disk, the CO line widths drop to values more 
consistent with the virial velocity of the host. In the most mas- 
sive 10'^* Mq, model Ql; Table 1) halo, this is manifested 
in broad {< a >^ 500-800 km s~') predicted CO emission 
line widths, whereas in the halos with lower circular velocity 
(quasars Q2, Q3), the mean CO line width is ~ 450 and 300 
km s"' , respectively. It is important to note that in all cases, 
while the aforementioned trends hold, a large range of line 
widths is permitted at all points as they are strongly sightline- 
dependent. Another way to view this is through the detailed 
distribution of line widths themselves. In Figure 9, we show 
a histogram of the sightline-dependent line widths during the 
quasar phase for models Q1-Q3. 

These derived linewidths are a natural consequence of our 
initial assumptions of quasar formation in massive halos, and 
that half the cold gas mass is in molecular phase. Within the 
confines of our initial assumptions, these results are consistent 
with virial arguments. Thus, for example, the broad emission 
lines seen in the massive ^ 10'^ M0 halo simply reflects the 
virial velocity of the host galaxy at the radial extent of the 
molecular gas distribution, which is of order ^ 550 km s"' . 

CO line measurements of J 1148-1-5251 by Bertoldi et al. 
(2003a), and Walter et al. (2003, 2004) showed CO (J=3- 
2, J=6-5 and J=7-6) emission lines with width ~ 280 ± 140 
km s"' (FWHM; corresponding to tr ~ 120 ± 60 km s"' for 
a Gaussian line). While this is a factor of 2-5 narrower 
than the median line width predicted in our models (Fig- 
ure 8 and 9), it is worthwhile to note the non-negligible 
fraction of sightlines for each of the quasar models that are 
compatible with the narrow observed line widths (Figure 9). 
Through the quasar phase and beginning of post-quasar phase, 

2-3% of sightlines in the most massive quasar host (Ql) 
have linewidths compatible with the line widths measured in 
Jl 148H-5251 (Bertoldi et al. 2003a; Walter et al. 2003, 2004). 
The smaller virial velocity of the lower mass halos (Q2 and 
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Fig. 5. — CO (J=3-2) emission from most massive quasar host (Ql) at z=6.73, while cold gas from the merger is still infalling. Top three panels are three 
orthogonal viewing angles, and the bottom six viewing angles are randomly drawn. The orientation for the line of sight is in the top right corner of each panel. 
Along many different viewing angles, multiple CO components are visible, similar to the observed morphology of Jl 148+525 1 (Walter et al. 2003), indicating a 
viable merger origin for Jl 148+525 1 . Each panel is four kpc on a side. The emission is in terms of velocity-integrated intensity (K-km s"' ), and the scale is at 
the bottom. 



Q3) naturally produce more sightlines with narrow line widths 
compatible with observations. During the quasar phase, the 
quasar formed in the 4 x 10'^ M0 halo (model Q2) shows 
CO (J=6-5) emission line widths consistent with observations 

5% of the time, and the quasar formed in the 1.5 x 10'^ 
Mq halo (model Q3) reveals CO (J=6-5) line widths consis- 
tent with observations 10% of the time. 

In an effort to understand the relationship of these models 
to the measured CO line width of Jl 148-1-5251, it is of inter- 
est to explore the origin behind the particular percentages of 
sightlines that are compatible with observations, and possible 



trends which may tend observations toward particular sight- 
lines. We thus focus the remainder of this section on this in- 
vestigation. We conclude the section with a discussion as to 
what circumstances may bring agreement between our simu- 
lations and observations, and the implications of these results 
with respect to potential observable tests motivated by these 
models. 

5.2. Ejfect of Merger Remnant Structure and Disk Formation 
on Line Widths 

We have thus far considered the formation of quasars hi- 
erarchically through multiple mergers whose physical condi- 
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Fig. 6. — Excitation dependent CO morpliologies: CO (J=l-0), (J=3-2) and (J=6-5) emission contours from most massive quasar liost (Ql) during the pealc of 
tlie quasar pliase. Columns are constant in redsliift, and rows are constant in transition. Multiple CO emission peaks are seen arising from merging cold clumps 
of gas, as well as circumnuclear gas in the (J=l-0) case. The gas in the central ~ kpc is highly excited, and thus does not show prominent (J=l-0) emission, but 
causes the (J=3-2) and (J=6-5) emission to be centrally concentrated. Multiple CO emission peaks are seen in higher lying transitions owing only to merging gas 
as the bulk of the gas in the central ~ kpc is highly excited. These results are robust across all modeled viewing angles. The viewing angle for each panel is 
the same, and is S = 0, </> = 0, for comparison with Figures 4 and 5. Each panel is four kpc on a side, and each row (transition) is on its own scale to facilitate 
interpretation. The scales are on the light, and the units are in terms of velocity integrated Rayleigh-Jeans temperature (K-km s"'). 



tions were derived self-consistendy from cosmological sim- 
ulations. It is possible that not all quasars at z 6 form via 
numerous violent merging events, but rather through a more 
'ordered' merger. It is thus worth quantifying the potential 
dependence of line width on merger history, and in particu- 
lar, how a more ordered merger remnant may affect the ob- 
served line widths. To provide a limiting example, we have 
conducted a test simulation of a coplanar binary merger in 
a halo of 10'^ Mq (model Q4; Table 1). The progenitors 
were initialized with a Hernquist (1990) profile, spin param- 
eter A=0.033 and circular velocity V200 ~ 600 km s~'. The 
disks had an initial gas fraction of 0.99, and the virial prop- 
erties were scaled to be appropriate for z=6 (Robertson et al. 
2006a). 

In Figure 10, we show the CO (J=3-2) morphology for 
the resultant quasar over two orthogonal viewing angles, and 
their corresponding unresolved emission spectra. The rem- 
nant forms a strong disk-like morphology, consistent with the 
findings of Robertson et al. (2006b) and Springel & Hern- 
quist (2005). As expected, the CO line width and profile is a 
sharp function of the viewing angle of the disk. Averaged over 



250 random sightlines, we find that the predicted CO (J=6-5) 
emission line width from this source is consistent with the ob- 
served line width J1148-I-5251 ^ 4% of the time, comparable 
to the more massive quasar hosts (Ql and Q2) which formed 
from multiple mergers. 

The percentage of sightlines compatible with observations 
in the binary merger Q4 arises from a limited range of angles 
a disk can be from face on to keep line widths within a partic- 
ular limit. Specifically, if one considers an inclined toy disk 
of pure gas with velocity dispersion o-yiriah then in order for 
observed velocity dispersions to fall below a particular value, 
(Jobs, the inclination angle from face on is limited by: 



< sm 



fobs 
f^virial 



(2) 



However, inclinations along both the polar (9) and azimuthal 
((p) axes have to be within this limit of a face-on configuration 
to keep the line of sight velocity dispersion below observed 
line widths. As such, the probability of having both 9 and (p 
randomly drawn such that they both fall below a critical value 
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Fig. 7. — CO (J=6-5) emission Unes from model quasars Ql, Q2 and Q3 at peak of their quasar phases viewed from three random sightUnes. The first row 
corresponds to quasar Ql, the second to Q2 and third row to Q3, The characteristic emission line widths drop with halo mass, though remain featured in each 
model with several substructure spikes superposed on broader Unes. Spectra have been convolved with a circiiar 5" Gaussian beam, and modeled at an angular 
diameter distance of 1 Gpc. 



to match an observed line width is 



PiO < Ooiit, (/> < (/icrit) = -T sin' 



1\ 2 



*^virial 



(3) 



where P « 4% for (Jobs = 180 km s"' (the upper Hmit a of 
the observed line in J1148+5251), and o"viriai = 600 km s~^ 
This probability then represents the upper limit of fraction of 
sightlines that will be compatible with the observed CO line 
width in J1148+5251 for disks with circular velocity ~ 500 
km s~' . 

This clarifies why the percentage of sightlines compati- 
ble with observations in the massive multiple-merger models 
(Q1-Q3) is relatively small. The cold gas in the three quasars 
which formed out of multiple non-idealized mergers (Q1-Q3) 
settles into rotating nuclear disks. As an aside, it is interest- 
ing to note that the amount of rotationally supported gas is 
seen to be dependent on the mass of the galaxy. Namely, the 
lower mass halos show a larger percentage of gas in stable ro- 
tation. To illustrate this, in Figure 11 we show the fraction of 
rotationally supported gas for quasars Q1-Q3 as a function of 
redshift, noting in particular the points of peak quasar activ- 
ity. When the most massive galaxy (Ql) is seen as a quasar, 
roughly 50% of the H2 gas is rotationally supported. Con- 
versely, in the lowest mass model (Q3), ^ 90% of the gas 
is rotationally supported. This may be a direct result of the 
amount of energy input from the central quasar during these 
times. Cox et al. (2007) demonstrated that the amount of 



rotationally supported gas in galaxy mergers decreases with 
increasingly efficient winds. During the peak of the quasar 
phase, the black hole luminosity is over an order of magni- 
tude brighter in model Ql than in Q3. 

In either case, though, large percentages of the gas in all 
of the quasars (Q1-Q3) are seen to be rotationally supported 
during the quasar phase. Because of this, the number of sight- 
lines compatible with observations is roughly characterized 
by Equation (3). Thus, the most massive model (Ql) shows 
line widths compatible with observations ~ 2-3% of the time, 
consistent with a predicted upper limit of ~ 4% for ~ 10^^ 
Mq halos. Similarly, in the lowest mass model (Q3), where 
~ 90% of the gas is rotationally supported, the 10% of sight- 
lines seen to be compatible with observed line widths is com- 
patible with the predicted upper limit of ^ 13%. The fact that 
the modeled fraction of sightlines is always slightly smaller 
than the toy model in Equation (3) owes to the fact that some 
of the gas in the galaxy is still highly dynamical and not viri- 
aUzed. 

5.3. CO Line Width-Quasar Luminosity Relation: Potential 

Selection Effects 

Given the dispersion of CO (J=6-5) line widths along dif- 
ferent Unes of sight, an interesting relation to explore is 
one between the optical quasar luminosity and CO emission 
line width in search of potential observational selection ef- 
fects which may tend observations of z ~ 6 quasars toward 
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Fig. 8. — Sightline-averaged velocity dispersion (tr) of CO (J=6-5) emis- 
sion lines (solid line, middle) as a function of redshift for quasar models Ql- 
Q3. The shaded area shows the range of derived cr over 250 randomly sam- 
pled sightlines, and the upper and lower dotted lines show the 2a linewidths 
for each snapshot. The spectra have been convolved with a circular 5" Gaus- 
sian beam, and binned to 50 km s"' . During the hierarchical buildup of the 
host galaxy, the H2 gas is highly dynamical, and the typical line widths are 
broader than the virial velocity of the host galaxy by a factor of ~ 1 .5 - 2. As 
the gas virializes during the quasar phase, the line widths drop, and roughly 
trace the virial velocity of the galaxy. A range of line widths are permitted 
throughout the evolution of the host galaxy, with larger numbers of sight- 
lines being compatible with the nan'ow (cr 120 km s"') detected line in 
J 1 148+5251 near the end of the quasar phase. The number of sightlines com- 
patible with observations naturally increases in the lower mass halos (Q2 and 
Q3) as the CO lines faithfully trace the virial velocity of the galaxy during 
the quasar phase. 



Fig. 9. — Histogram of sightline-dependent line widths during quasar phase 
for three halo mass models. The mean line width traces the virial velocity of 
the host galaxy, and thus become narrower in lower mass galaxies. In all 
cases, a broad range in linewidths is observable. 
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Fig. 10. — CO (J=l-0) morphology and unresolved line profiles of rem- 
nant formed from binary merger simulation (model Q4) across two orthog- 
onal sightlines (face-on and edge-on). The progenitors were initialized on a 
coplanar orbit in a ~ 10'^^ M0 halo. The resultant quasar host galaxy has 
large rotating molecular disk, resulting in narrow observed emission lines in 
face-on viewing angles, and broad lines in edge-on viewing angles. Approxi- 
mately 4-5% of sightlines show narrow line widths compatible with observa- 
tions of J1148+5251 (see Equation 3). The maps are 12 kpc across, and the 
scale is in units of K-km s"' . 

smaller line widths. In Figure 12, we plot the a from the 
CO (J=6-5) line width as a function of the attenuated rest- 
frame B-band luminosity over 5000 lines of sight throughout 
the quasar phase for the most massive model, Ql. We utilize 
the methodology of Hopkins et al. (2005a) in computing the 
dust-attenuated quasar luminosity, and include contributions 
from both the stellar component, as well as the central AGN. 

There is a general trend for sightlines which show the 
brightest rest-frame Z?-band luminosity to have smaller CO 
line widths. This can be understood via decomposition of 
the quasar luminosity into its stellar and AGN components. 
Specifically, while the stellar luminosity does not vary much 
with viewing angle, the contribution of the central black hole 
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Fig. 11. — Fraction of rotationally supported molecular gas in quasars Ql- 
Q3. Gas is considered to be in rotation^ support if its rotational velocity is 
at least 80% the expected circular velocity at that radius. Arrows point to 
redshift of peak quasar activity (also listed in Table 1). Gas is seen to more 
easily settle into a stable rotational configuration in the lower mass quasar 
hosts, though all quasars naturally form relatively strong disks during the 
quasar phase. 

to the total luminosity is strongly dependent on the viewing 
angle with respect to the rotating molecular gas. Directions 
which view the molecular disk in an edge-on configuration 
result in a heavily obscured central AGN. Along these sight- 
lines, the emergent CO emission line is typically broad (e.g. 
bottom right panel. Figure 10). Conversely, when the molec- 
ular disk is seen in a more face-on viewing angle (and thus 
has narrower CO emission lines), the central AGN can be 
viewed relatively unobscured, and the attenuated rest-frame 
B-band luminosity is consequently higher. This effect is typi- 
cal during the quasar phase, and only rarely is the black hole 
relatively unattenuated through an edge-on sightline, which 
causes broad CO lines to be visible when the rest-frame B- 
band luminosity peaks. 

The results of this relationship between bolometric lumi- 
nosity and line width suggest a potential selection effect 
which may cause quasars selected for optical luminosity to 
have systematically lower CO line widths, owing to the pre- 
ferred face-on viewing angle for the molecular disk. In Fig- 
ure 13, we plot the percentage of sightlines with line widths 
compatible with observations as a function of limiting rest- 
frame B-band luminosity. For the highest flux cuts, the frac- 
tion of sightlines with narrow line widths increases from 2-3% 
to ~ 10%. This selection effect is robust across both the lower 
mass halo models as well. In models Q2 and Q3, the fraction 
of sightlines compatible with observations increases to ~ 15 
and 25% respectively, for the highest luminosity cuts. 

Finally, we note that the viewing angles corresponding to 
the smallest (2a) line widths in our simulations typically fall 
within a 25° range in polar and azimuthal angle from face-on. 
This range of angles is consistent with recent observational 
studies which have suggested that the commonly observed 
narrow CO line widths of high-z quasars may correspond to 
a preferred viewing angle of 10 — 15° from face-on (CarilU & 
Wang 2006; Wu 2007). 

5.4. Interpretation of z ^ 6 Quasar Observations 

Our simulations show that z ~ 6 quasars which form in 
massive halos will characteristically have broad mean line 
widths, consistent with simple virial arguments. The mean 
simulated fine widths are in apparent contradiction to ob- 
served narrow (a ^ 120 ± 60 km s~^) line widths of the sole 
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Fig. 12. — CO (J=6-5) emission line width (it) versus attenuated rest-frame 
B-band luminosity. The brightest quasar luminosities correspond to face-on 
molecular disk configurations where the central AGN is the least obscured. 
These face-on sightlines additionally show nanow CO line widths. Thus, 
surveys which select quasars for optical luminosity may preferentially select 
objects that have narrow CO line widths. This may have important conse- 
quences for using CO as a dynamical mass indicator in quasars. 




^ 10" 10'^ 
Limiting B-band Luminosity (attenuated) 

Fig. 13. — Percent of sightlines of ~ 10'^ Mq quasar over which line 
widths are compatible with observations versus limiting rest-frame B-band 
luminosity. Natural molecular disk formation in the quasar host galaxy gives 
rise to a selection effect which enhances the likelihood that flux-limited op- 
tical surveys will view the quasar in a face-on configuration. These viewing 
angles also typically have narrower CO line widths, and thus may increase the 
chances of viewing a narrow emission line in quasars which is not necessarily 
characteristic of the true sightline-averaged mean line width. Line widths are 
derived from CO (J=6-5) line profiles, and compared to upper limit of sole 
CO detection at z ~ 6, Jl 148+525 1 . Rest-frame B-band luminosities are dust 
attenuated, and calculated using the methodology of Hopkins et al. (2005a). 

CO detection at z ~ 6 (Walter et al. 2004), though a non- 
negligible fraction of sightlines (ranging from 2-3% for the 
most massive model to ~ 10% in the lowest mass model) are 
compatible with this observation. Potential selection effects 
owing to molecular disk formation in these galaxies will in- 
crease the probability of narrow-line detection when select- 
ing quasars for optical luminosity, ranging from ^ 10% in the 
most massive quasar host to 25% in less massive ones. The 
full dispersion in line widths is predicted to become more ap- 
parent at lower optical luminosities (Figure 12). 

Our models find that quasars which form in the lower end 
of our halo mass range at z ~ 6 (e.g. quasars Q2 and Q3) 
may have similar rates of detection as those which form in 
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the most massive ^ lO'"' M0 halos (e.g. quasar Ql). This 
owes to the competing effects of smaller quasar lifetimes in 
the lower mass halos, but more lower mass halos in the sim- 
ulated cosmological volume. Specifically, the quasars formed 
in the 4 X 10^^(1.5 x 10^^) Mq halos have luminosities > lO'^ 
L0 for a factor of ~ 2.5(14) less time than the quasar formed 
in the 10'^^ Mq halo (Table 1). However, standard halo mass 
functions predict more low mass halos than massive ^ lO'"' 
Mq halos (Press & Schechter 1974; Sheth & Tormen 2002; 
Springel et al. 2005). The cosmological simulations of Li 
et al. (2007) (which were used in this work) found approxi- 
mately 3(7) times as many halos of mass 5 x 10'^(2 x 10'^) 
Mq, compared to the single ~ 10'-^ Mq halo identified in the 
simulation box at z 6 (Li et al. 2007, Figure 14). In this 
sense, quasars which form in lower mass halos may be identi- 
fied at similar rates as those which form in ~ lO'^ Mq halos. 

One result of this work is to motivate observational tests of 
these models. A direct prediction of these simulations is that 
a CO survey from a large sample of quasars at z ^ 6 which 
probes lower on the optical luminosity function may directly 
constrain the range of potential line widths which originate 
from high redshift quasars. At the median rest-frame Z?-band 
luminosity, our models predict that a large range of line widths 
should be observed (Figure 12). A potential caution associ- 
ated with this test is that physical processes on scales below 
the resolution of our simulations may limit optical detections 
of quasars with edge-on disks at z ^ 6. For example, if a 
dusty molecular torus exists on scales smaller than ^ 100 pc 
and provides high levels of obscuration along sightlines other 
than face-on, then the variation of optical luminosity with in- 
cUnation angle will be steeper than suggested by Figure 12. It 
is thus not a straightforward assumption that optical surveys 
will need to probe only an order of magnitude lower in rest- 
frame B-band luminosity to test these models as Figure 12 
suggests, but rather it is in the limit that quasars with rela- 
tively incUned disks at z ~ 6 can be detected that these mod- 
els predict a broad range of CO line widths at lower optical 
luminosities. 

A more clear test may come from observations of the pro- 
genitors of z ~ 6 quasars themselves. For example, a direct 
prediction from these models is that CO observations of ei- 
ther the most massive progenitor galaxies prior to the merger 
(e.g. at z ^ 8; see Table 1 of Li et al. 2007), or the ongo- 
ing mergers themselves may exhibit a large dispersion in CO 
line widths, with median velocity dispersion reflective of the 
halo virial velocity (Figure 8). TTie identification of potential 
progenitors at z >7 is predicted to be feasible through z-band 
dropouts (Robertson et al. 2007). 

In either case, these models suggest that surveys at z ^ 6 
which observe either the most massive progenitors of quasar 
host galaxies, or highly incUned disks (closer to edge-on, and 
likely lower luminosity) associated with z ^ 6 quasars will see 
a broader dispersion of CO line widths. The model halos pre- 
sented here reflect the range of halo masses in our cosmo- 
logical simulation which were feasibly able to create a z ~ 6 
quasar. As such, based on the line widths seen in our low- 
est mass host galaxy (Q3), samples of CO detections at z ~ 6 
which probe quasars with a range of disk inclination angles 
should find a median line width at least ^ twice the value of 
the sole detection. Surveys which observe quasars with in- 
cUned molecular disks and still find consistently narrow line 
widths may reflect an inability of our radiative transfer simu- 
lations to fully capture the appropriate physics necessary to 
predict accurate CO fine widths form these early Universe 



galaxies. If, for example, our assumptions regarding molec- 
ular gas content or CO abundances are incorrect, it could be 
that CO emission is preferentially seen in lower velocity gas 
in the host galaxy. Alternative possibilities may include super- 
Eddington accretion for the central black hole which may al- 
low for luminous quasars at z ~ 6 in less massive halos than 
those adopted here (e.g. Volonteri & Rees 2005). It is, how- 
ever, attractive that the simulations presented here do provide 
a model for the CO emission from the earliest quasars which 
are consistent with observations of the CO excitation and mor- 
phology of Jl 148-1-5251. More observations to fully deter- 
mine the nature of CO line widths in z ^ 6 quasars will be 
necessary to assess the vahdity of this aspect of our model- 
ing. 

6. SPHEROID, BLACK HOLE, AND DYNAMICAL MASS 

Because submillimeter-wave radiation typically does not 
suffer the heavy extinction characteristic of optical emission, 
CO lines are often used as dynamical mass indicators in dusty 
starburst galaxies. Dynamical masses are derived assuming 
that the emitting gas is rotationally supported, and that the 
line width provides a measure of the rotational velocity. Here 
we assess the usage of CO-derived dynamical masses in z ~ 6 
quasars. As a case study, we will focus on the most mas- 
sive simulation (^^ 10^^ Mq; model Ql), though the trends 
are generic for all quasars presented in this work. 

Generally, theoretical arguments have predicted that the 
Mbh-<^v and MfiH-A/buige relations show only weak evolution 
with redshift (Robertson et al. 2006a). Using numerical sim- 
ulations of galaxy mergers by Cox et al. (2006b) and Robert- 
son et al. (2006a) that include black hole feedback, Hopkins 
et al. (2007a) found that the normalization of the MfiH-A^buige 
relation shows weak (~ 0.3 - 0.5 dex) trends toward larger 
MfiH/Mbuige from z= - 6. The simulations presented here 
provide additional support for a scenario in which the stel- 
lar bulge and central black hole grow coevally in the earliest 
galaxies. 

In the most massive galaxy presented here, the supermas- 
sive black hole grows rapidly during the hierarchical buildup 
of the quasar host galaxy, and reaches a total mass of ~ 
2 X 10^ M0 during the peak quasar phase, similar to black 
hole mass estimates in J1148H-5251 (Willott et al. 2003). Ow- 
ing to the extreme star formation rates which can be as large as 
^ 10"* M0yr~' between redshifts z=9 and z=8 during the final 
violent mergers (e.g. Figure 3), the bulge reaches a total mass 
of ~ 10^^ M0 by the quasar phase. The black hole and stel- 
lar bulge masses are related such that Mbh = 0.002 x Mbuige 
during the quasar phase, roughly consistent with the local 
MsH-Mhuige relation (Li et al 2007; Magorrian et al. 1998; 
Marconi & Hunt 2003). We note, though, that during the peak 
of the quasar activity, the Mbh - fv relation is not necessarily 
obeyed as the stellar bulge is not dynamically relaxed. 

The massive starburst at z >7 results in an ISM highly en- 
riched with metals, consistent with the observed Fell/Mgll 
abundances, [CII] emission CO emission, and dust masses 
for z ~ 6 quasars (Barth et al. 2003; Bertoldi et al. 2003b; 
Dietrich et al. 2003a;, Freudling et al. 2003; Jiang et al. 
2006; Maiolino et al. 2005; Walter et al. 2003, 2004). Li 
et al. (2007) have found that the metallicity in the simulated 
~ 10^^ Mq quasar presented here is solar to supersolar dur- 
ing the quasar phase, owing to the ~ lO'' M0yr~^ starburst 
during the hierarchical merging process of the quasar forma- 
tion. These findings are consistent with the mean abundances 
of ~ 4 times solar in a sample of 4 <z< 5 quasars observed by 
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Dietrich et al. (2003b) and imply that black hole growth and 
stellar bulge formation are correlated at high redshifts. These 
simulations support findings that the relationship is regulated 
by feedback from supermassive black holes (e.g. Di Matteo, 
Springel & Hernquist 2005; Robertson et al. 2006a; Hopkins 
et al. 2007a). 

CO observations of high-z quasars have suggested that 
central black hole masses may be excessively large com- 
pared to the stellar bulge mass as predicted by the present 
day MBH-A^buige relation. Dynamical mass estimates of 
J 1 148+5251 using CO line measurements have indicated that 
the stellar bulge may be undermassive by a factor of 10-50 
if the present day Msn-Mbuige relation holds at z ^ 6 (Walter 
et al. 2004). Studies of other quasars at z ^3 using CO line 
widths as a proxy for enclosed mass have arrived at similar 
conclusions (Shields et al. 2006). 

In order to investigate the usage of CO as a dynamical mass 
indicator, in Figure 14, we plot the evolution of the median 
dynamical mass of the quasar host galaxy from the most mas- 
sive simulation (Ql) derived from the CO (J=6-5) emission 
line over 250 random lines of sight, as well as 3(7 contours for 
the observed range of linewidths. In particular, we show 

Md,„ = k—-— (4) 
G sm(i) 

where we use an adopted k of 8/3 and assume an average incli- 
nation angle of 30°. We plot the total mass of the host galaxy 
(within the central 2 kpc) throughout its evolution, including 
the multi-phase ISM, dark matter, stars, and black holes. We 
additionally plot the ratio of the median derived dynamical 
mass to the total mass enclosed. The dynamical mass is de- 
rived from the width of the CO (J=6-5) line as the CO excita- 
tion peaks at the J=6 level, and this transition thus best trace 
the bulk of the molecular gas during the quasar phase of our 
simulated galaxy. We note, however, that the line widths from 
lower transition lines do not deviate much from these trends. 

During the hierarchical buildup of the host galaxy (pre- 
quasar phase), the gas is still highly dynamical, and not com- 
pletely virialized. This results in large CO line widths, and 
consequently derived dynamical masses which largely over- 
estimate the true mass. During this time, the typical derived 
dynamical mass results in an overestimate of the true mass 
by a factor of ~ 2-5. As the gas becomes more rotationally 
supported during the quasar phase, the CO (J=6-5) line width 
serves as a better tracer of the mass enclosed. These results 
are consistent with similar derivations by Greve & Sommer- 
Larsen (2006) who found that CO can serve as an accurate 
dynamical mass tracer in merger simulations to within 20%. 

An important point from Figure 14 is that a large sightline- 
dependent range of values are possible for derived dynami- 
cal masses. While the median line widths provide reason- 
able estimates of the true mass, many lines of sight permit 
significant underestimates. The radiative transfer simulations 
presented here may then bring some resolution between CO 
observations which are suggestive of a strong evolution in 
the MBH-A^buige relation at high redshift, and models which 
imply a lesser evolution between present epochs and early 
times. Figure 14 shows that a large range of dynamical masses 
may be inferred simply based on observed viewing angle. 
While the median CO line widths accurately traces the en- 
closed mass for the bulk of the quasar phase, potential se- 
lection effects (Figures 12 and 13) may bias masses derived 
from CO line widths from quasars toward an underestimate 
of the true mass. These model results are consistent with the 





Fig. 14. — Top: Median dynamical masses derived from line widths 
over 250 random sightlines through the hierarchical buildup and evolution 
of quasar host galaxy (red crosses). The shaded region represents the 3cr 
range of sightline-dependent derived dynamical masses, and the solid line 
the true total mass within the central 2 kpc. The dynamical masses are de- 
rived assuming an incHnation angle of 30° , and a molecular spatial extent of 
2 kpc. Bottom: Ratio of median derived dynamical mass to total mass. Dur- 
ing the hierarchical buildup and early quasar phase, much of the gas is highly 
dynamical; consequently, derived dynamical masses from CO line widths do 
a poor job representing the true enclosed mass, typically overestimating by 
factors of 2-5. As the gas virializes the sightline-averaged CO line width is a 
better estimator of the dynamical mass. A large range of values is permitted, 
consistent with the spread in CO line widths seen in Figure 8. 

recent study by Wu (2007) who found that CO line widths 
provided a particularly poor estimate of bulge velocity disper- 
sion for local Seyferts with line widths narrower than ~ 400 
km s"' (FWHM; cr ^ 170 km s"'). In contrast, the inclination- 
corrected CO line width was found to correspond well with 
the bulge velocity dispersion. 

7. COMPARISONS TO OTHER HIGH REDSHIFT POPULATIONS 

A natural question which arises from this work is, how do 
the CO properties of these extreme 6 objects relate to 
other known starburst and AGN populations? Two extreme 
high redshift galaxy populations which may serve as interest- 
ing comparisons are: 1. Optically selected quasars; and 2. 
Dusty submillimeter selected galaxies at z ~ 2 (see respective 
reviews by Solomon & Vanden Bout, 2005, and Blain et al. 
2002). 

7.1. Quasars 

Only a handful of quasars have been detected in molecular 
line emission owing to beam dilution and limited sensitivity at 
millimeter and submillimeter wavelengths. As of the writing 
of the recent review by Solomon & Vanden Bout (2005), there 
have been about 16 quasars at z >1 for which molecular gas 
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emission properties have been pubhshed. 

Deciphering the molecular gas morphology remains an is- 
sue for most quasars since the majority of them have been 
detected in molecular line emission with the help of gravi- 
tational lensing, resulting in multiple imaging. The spatial 
extent of the molecular emission in most quasars appears to 
range from ^1-5 kpc (Solomon & Vanden Bout, 2005), sim- 
ilar to the simulations presented here as well as observations 
of Jl 148H-525 1. Multiple CO emission peaks are detected less 
often, though it is not clear whether or not this is a spatial res- 
olution issue. Of seven imaged quasars listed in the recent 
review by Solomon & Vanden Bout (2005), three have clear 
CO companions (including J1148H-5251), with distances in 
the emission peaks ranging from 1.7-8.7 kpc (e.g. Carilli et 
al. 2002). The more sensitive, and higher spatial resolution 
observations that will be routinely achieved with ALMA will 
clarify the CO morphology of quasars. 

The CO excitation conditions that have been measured in 
quasars have primarily been at z ^ 2. Barvainis et al. (1997) 
find that the CO is substantially excited in the Cloverleaf 
quasar (H1413-I-117, z~ 2.5) through the 1=1 level. High 
spatial resolution observations of two z ^ 4 quasars by Car- 
illi et al. (2002) find similar excitation conditions in CO, both 
of which have multiple CO emission peaks in their morphol- 
ogy. More recent large velocity gradient modehng of z ~ 4.7 
quasar BR 1202-0725 also found a CO SED that peaks at J=7 
(Riechers et al. 2006b). The high excitation conditions ob- 
served in quasars are indicative of extremely warm and dense 
molecular gas heated by ongoing star formation, and are con- 
sistent with active star formation during the quasar phase. 

The quasar sample of Solomon & Vanden Bout reports 
CO line widths in high-z quasars ranging from roughly a ^ 
100-250 km s"' (Carilh & Wang 2006). Large bandwidth 
observations with the 4GHz COBRA correlator on OVRO by 
Mainline et al. (2004) of a sample of three z=2-3 quasars 
showed similarly narrow CO line width. 

While these are in apparent contrast to the large mean 
line widths in z ^ 6 quasars reported here, the narrower line 
widths of quasars from 2 < z < 5 may be simply explained by 
the evolution of Vc with decreasing redshift. Namely, while 
quasars from 2 < z < 5 appear to form in halos of compara- 
ble mass to those at z 6 (Croom et al. 2005; Shen et al. 
2007), the expected circular velocity will naturally be lower 
than their z ~ 6 analogs by a factor of y'(l+z)/(l+6). As 
a result, the circular velocity for these halos at z ~ 2 would 
be of order a ~ 300-350 kms"\ consistent with cluster- 
ing measurements of z ~ 2 quasars (e.g. Croom et al. 2005; 
Porciani & Norberg 2006; Myers et al. 2006; Hopkins et 
al. 2007d) and studies of the quasar proximity effect (e.g. 
Faucher-Giguere et al. 2007; Kim & Croft 2007; Nascimento 
Guimaraes et al. 2007) , and slightly larger than measured 
CO line widths (Solomon & Vanden Bout, 2005), though it is 
possible that there is a potential selection bias toward quasars 
with narrow fine widths (§ 5.3) at lower redshifts as well. 

In this context, we emphasize that the simulations presented 
here do not predict large (cr ~ 500-800 km s~^) line widths 
for quasars which form at lower (e.g. z 2) redshift, even if 
they form in halos of comparable mass to luminous sources 
at z 6. To illustrate this, we examine the predicted CO line 
widths of the binary merger simulation of Narayanan et al. 
(2006a; also presented in Robertson et al. 2006a). The cir- 
cular velocity of the progenitor galaxies was Vc=160 km s"'. 
The total masses of the progenitor galaxies in this example 



was 4.8 X lO" h"'M0, and the final merger produced a central 
black hole mass of ^ 5 x 10^ h"'M0. The a of the unresolved 
CO (J=l-0) fine profile ranges from 100-150 km s"' during 
the peak of the quasar phase, over three orthogonal viewing 
angles, consistent with CO fine width measurements of z ~ 2 
quasars, and reflective of the host virial velocity. 

These models then suggest a self-consistent picture which 
naturally explains the evolution of CO line widths as function 
of redshift in terms of the circular velocity of the quasar host 
galaxy halo. Future modehng will need to quantify the ex- 
tent of potential molecular disk-driven selection effects (e.g. 
§ 5.3) for quasars at lower redshifts (z ~ 2-4). Indeed other 
works have suggested that quasars in this redshift range may 
be observed at a preferred shallow range of viewing angles as 
evidenced by their CO line widths (e.g. CariUi & Wang 2006). 

7.2. SubmilUmeter Galaxies 

Submillimeter galaxies (SMGs) represent a class of mas- 
sive, dusty starbursts at z ~ 2 largely detected by blind sur- 
veys with the SCUBA and MAMBO bolometers. Typical in- 
frared luminosities in these sources of ~ 10^' L0 correspond 
to a SFR of > 1000 MQyr"', assuming an insignificant AGN 
contribution (Small 2006). X-ray and IR studies have shown 
that these galaxies are known to contain embedded AGN, al- 
though their relative contribution to the bolometric luminosity 
is uncertain (Alexander et al. 2005a,b; Donley et al. 2005; 
Polletta et al. 2006). Optical morphologies of SMGs indi- 
cate that many are interacting and/or mergers (Chapman et al. 
2004). Recent CO morphologies and emission line profiles 
have furthered this scenario (Greve et al. 2005; Tacconi et 
al. 2006). Many studies have pointed to a picture in which 
SMGs may be ongoing mergers at z ^ 2, though in a pre- 
quasar phase (e.g. Blain et al. 2002). This combined with the 
fact that SMGs are the most massive and actively star forming 
galaxy population at z ~ 2 make SMGs an interesting compar- 
ative for our simulated z ~ 6 quasar. 

The excitation characteristics have been observed in multi- 
ple CO emission lines for only one case: z=2.5 SMG SMM 
J16359H-6612. In this galaxy, the CO flux density peaks at J=5 
(WeiB et al. 2005a), consistent with highly excited gas. The 
excitation conditions are similar to those seen in our simula- 
tions when the SFRs are comparable. 

The molecular morphology in SMGs closely resembles the 
extended CO emission seen in simulations and observations 
of Jl 148H-5251, as well as in the pre-quasar phase galaxy in 
our simulations at z=7-10. The average CO FWHM radius in 
the recent Tacconi et al. (2006) sample is 2 kpc. The large 
spatial extent of CO emission in SMGs combined with an 
apparent lack of strong quasar activity has been interpreted 
as being the consequence of extremely massive and gas rich 
galaxy mergers early in the evolution of the galaxy, and prior 
to the quasar phase (e.g. Tacconi et al. 2006). Indeed, if 
SMGs form through hierarchical mergers, then the models 
presented here suggest that the spatially extended and dis- 
turbed CO morphologies seen in SMGs further indicate that 
these galaxies may be mergers prior to an optical quasar phase 
(e.g. Figure 4). 

The CO line widths of SMGs are more enigmatic. The aver- 
age CO emission line from an SMG is typically about double 
the width of that from a quasar of similar redshifts (e.g. Carilh 
& Wang 2006; Greve et al. 2005; Tacconi et al. 2006). The 
average line width in SMGs is cr - 330 km s'^ (FWHM ~ 780 
km s"^), compared to a mean cr of 130 km s"' seen in z ~ 2 
quasars (Greve et al. 2005). Within the context of the models 
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presented in this work, a number of physical motivations for 
these differing line widths may be at play. 

First, in our models the mean Une widths from galaxies with 
virialized cold gas are seen to roughly correspond with the 
circular velocity of a host halo. When comparing galaxies 
of similar redshifts, if SMGs are dynamically relaxed, this 
would imply that SMGs may originate in more massive ha- 
los than quasars. Indeed, clustering measurements made by 
Blain et al. (2004) have suggested that SMGs are hosted by 
halos slightly larger than quasars at comparable redshifts. If 
SMGs reside in halos >4 times more massive than typical 
quasar host galaxies, the difference in CO Une widths may be 
accounted for. 

An alternative explanation for the line widths may arise 
from an evolutionary standpoint. Specifically, as Figure 8 
shows, there is a typical drop in CO line width by a factor 
of ^ 2 close to the quasar phase of a merging galaxy sys- 
tem. This owes to gas in the violent environs of a galaxy 
merger becoming rotationally supported late in the evolution 
of the merger. If SMGs are a class of objects hosted by ha- 
los of similar mass to their quasar counterparts, then it may 
be that SMGs are simply massive mergers at z 2 prior to 
their active quasar phase. Certainly CO line profiles and 
morphologies from SMGs are consistent with merging ac- 
tivity (Greve et al. 2005; Tacconi et al. 2006; Narayanan 
et al. 2006; D. Narayanan et al., in prep). Moreover, nu- 
merical simulations of merging galaxies coupled with self- 
consistent radiative transfer solutions have pointed to a pic- 
ture in which SMGs may be mergers caught in a phase of 
massive black hole growth, though prior to an optical quasar 
phase (Chakrabarti et al. 2007b). 

Other authors have suggested that SMGs may be similar to 
quasars in their place in galaxy evolution, but simply viewed 
more edge-on. Carilli & Wang (2004) suggested that if these 
2-^2 galaxy populations are truly of the same class, then one 
can infer a mean viewing angle of ~ 13° for quasars, whereas 
SMGs are more likely randomly oriented. 

Finally, it may be that the CO emission properties of SMGs 
are not explained by physical models such at those presented 
in this work. If, for example, the massive star formation rates 
are fueled at least in part by accretion of gas from the host 
halo (which is seen to occur for at least some galaxies in cos- 
mological simulations; see e.g. Finlator et al. 2006), the CO 
line widths may not reflect the circular velocity of the system, 
at least during phases of elevated star formation, as they do in 
merging galaxies. 

In summary, correspondence between some observed 
trends of molecular line emission in the models of hierar- 
chical z ^ 6 quasar formation and SMGs which form at later 
times suggest that it is plausible that SMGs fall naturally into 
a merger-driven evolutionary sequence, though, at least from 
molecular line diagnostics alone, its location on this sequence 
is not completely clear. Of course, alternative scenarios can- 
not be ruled out here. Further models of hierarchical galaxy 
mergers appropriate for z ~ 2 will have to be examined in or- 
der to further quantify the relation of CO emission properties 
to the evolution of SMGs. 

8. SUMMARY AND CONCLUSIONS 

We have apphed non-LTE radiative transfer simulations 
to cosmological and hydrodynamic galaxy formation simu- 
lations to predict the CO emission from representative z ^ 6 
quasars that are modeled to form hierarchically in massive 
lO^'^- 10^^ Mq halos. We made predictions concerning the 



CO excitation patterns, morphologies, and line widths in this 
extreme class of objects. We further made broad compar- 
isons to the only current CO detection at z >6, Jl 148H-525 1 at 
z=6.42. Our main results are the following: 

1 . Owing to very warm and dense conditions in the molec- 
ular ISM, the CO flux density is predicted to peak at the 
J=8 level during the early, hierarchical formation pro- 
cess of the quasar host galaxy (z ^ 8), when the SFR 
can be as high as ^ 10"* M0yr~^ During the peak 
quasar phase, the central AGN reduces the nuclear star- 
burst, and the SFR drops to ^ 10^ Moyr"'. Conse- 
quently, the CO flux density peaks at the J=5-6 level. 
These excitation conditions are indicative of an ongo- 
ing starburst, and are consistent with observations of 
J1148-H5251. As the gas becomes more diffuse and the 
starburst dies down in the post-quasar phase (z <6), the 
peak in the CO flux is predicted to drop to Jw 3. 

2. The CO morphologies of z ~ 6 quasars may exhibit 
multiple emission peaks during the active quasar phase, 
owing to separated peaks of high density emission that 
have not yet coalesced. The multiple emission peaks in 
the morphology of the CO (J=3-2) gas during the quasar 
phase of our simulations are very similar to observa- 
tions of Jl 148H-5251, and are robust along many view- 
ing angles. These results imply that a merger-driven 
formation scenario for z ~ 6 quasars produces CO mor- 
phologies consistent with that of Jl 148-H525 1. 

3. On average, the CO line widths from quasars at z ^ 6 
are reflective of the virial velocity of the host halo, 
though there exists a large sightline-dependent disper- 
sion in Une widths. During the hierarchical buildup 
of the host galaxy, the median line widths are roughly 
twice the virial velocity, and settle to the virial velocity 
during the quasar phase. During the quasar phase, the 
sightline- averaged line width for the ~ lO'^ M0 halo is 
~ 500-800 km s'K In the lowest mass halo (~ 10'^ 
Mq) the sightline averaged line width is ~ 300 km s"^. 

4. A fraction of sightlines in each model is compatible 
with observations, and is a strong function of halo mass. 
Specifically, the number of sightUnes with narrow line 
widths compatible with observations increases with de- 
creasing halo mass. The most massive 10^^ M0 halo 
shows ^ 2-3% of sightlines compatible with observa- 
tions. The ^ 10'^ Mq halo has line widths similar to 
observations of J1148H-5251 ~ 10% of the time. The 
percentage of sightlines compatible with observations 
increase owing to selection effects (next point). 

5. Quasars at z 6 selected for optical luminosity may 
preferentially be in a face-on configuration as this pro- 
vides the least obscuration of the central black hole. In 
these configurations, the CO Une widths are narrower, 
thus causing quasars selected for optical luminosity 
to preferentially have narrower line widths than their 
sightline averaged values. The fraction of sightlines 
with line widths compatible with observations increases 
to 10-25% when considering quasars selected for opti- 
cal luminosity. This suggests that these models may 
be in agreement with observations if Jl 148H-5251 is be- 
ing observed in a face-on configuration. A direct con- 
sequence of these selection effects is that in order to 
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observe the full dispersion in CO line widths in z ^ 6 
quasars, observations must probe quasars with edge-on 
molecular disks (which are typically lower on the opti- 
cal luminosity function). 

6. Because of the evolution of Vc with redshift of halos 
of similar mass, quasars which form in ^ lO'^- lO'^ 
M0 halos (as they are thought to; see e.g. Croom el al. 
2005; Shen et al. 2007) at lower redshift will naturally 
have smaller line widths, compatible with observations. 
We explicitly show this by examining the line widths 
of a binary merger simulation appropriate for present 
epochs. In this light, the simulations presented here do 
«of predict large (cr ^ 500-800 km s"') CO line widths 
for quasars which form at lower redshifts, but rather a 
suggestive self-consistent model for the potential ori- 
gin of CO line widths in quasars at both low and high 
redshift. 

7. Our merger-driven model for quasar formation predicts 
a host galaxy that lies on the MBH-A^buige relation during 
the active quasar phase (Li et al. 2007; Robertson et al 
2006a; Hopkins et al. 2007a). During the hierarchical 
buildup of the host galaxy, the median CO line width 
tends to typically overestimate the dynamical mass by 
a factor of 2-5 as much of the gas is highly dynamical, 
and not viriahzed. During the quasar phase, dynamical 
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masses derived from the median line widths are a better 
representation of the true mass. There is a large range 
in derived dynamical masses coincident with the large 
sightline-dependent range of line widths seen at a given 
time. If selection effects are in place such that molecu- 
lar disks in observed high-z quasars are typically close 
to face-on, CO-derived dynamical masses will prefer- 
entially underestimate the true mass unless the shallow 
viewing angle is accounted for. 
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